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The Sun is a source of high energy neutrinos (E > 10 GeV) produced by cosmic ray interactions 
in the solar atmosphere. We study the impact of three-flavor oscillations (in vacuum and in matter) 
on solar atmosphere neutrinos, and calculate their observable fluxes at Earth, as well as their event 
rates in a kilometer-scale detector in water or ice. We find that peculiar three-flavor oscillation 
effects in matter, which can occur in the energy range probed by solar atmosphere neutrinos, are 
significantly suppressed by averaging over the production region and over the neutrino and antineu- 
trino components. In particular, we find that the relation between the neutrino fluxes at the Sun 
and at the Earth can be approximately expressed in terms of phase-averaged "vacuum" oscillations, 
dominated by a single mixing parameter (the angle #23). 
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I. INTRODUCTION 

The Sun is a well-known source of low-energy electron 
neutrinos (E < 20 MeV) , continually produced in its core 
by nuclear reactions. The study of this flux has been one 
of the most exciting enterprises in astroparticlc physics, 
and has provided us with crucial evidence for v flavor 
oscillations (see, e.g., 0,0 and references therein). 

High-energy neutrinos and antineutrinos (with E > 
10 GeV) are also continually produced by cosmic ray in- 
teractions in the solar atmosphere [1, 3, [H, [H, 0] , which 
lead to the production of secondary mesons and their 
subsequent decay into v e , v e , 1/^, and collectively 
referred to as "solar atmosphere neutrinos" (SAz^) here- 
after. From pion production dominance, at the source the 
SAis flavor ratio is expected to be (f) e : <\>^ : <fi T ~ 1 : 2 : 0, 
with approximately equal fluxes of v and V. Concerning 
the SAu energy spectrum, the Monte Carlo simulation in 
Q finds a power-law decrease (~ E~ 3 ) for E > 100 GeV. 
At lower energies, heliomagnetic effects are expected to 
make the spectrum less steep (as ~ E~ 2 ) @, @]. 

At the Earth, within the solid angle subtended by the 
Sun, the SAi^ flux significantly exceeds the more familiar 
"atmospheric neutrino" flux generated by cosmic rays in- 
teracting with the Earth atmosphere. In fact, because of 
the different scale heights, the density at the first inter- 
action point is lower in the solar atmosphere than in the 
terrestrial one, and therefore a larger fraction of mesons 
can decay into neutrinos instead of being absorbed. The 
expected SAv signal can be quantified in terms of O(10) 
events per year in a km 3 of water (or ice) above a thresh- 
old E > 100 GeV @, and might thus be detected by 
large- volume v telescopes like IceCube ||. 

The characterization of the SAz/ signal is an impor- 
tant goal, both in itself and because it might repre- 
sent a "background" in searches for signatures of pos- 
sible weakly interacting massive particle (WIMP) anni- 
hilation in the Sun @, |lo|, El EI El El, El ■ (Predicted 



WIMP event rates are rather model -dep endent , ranging 
from 10~ 12 to 10 3 events/year/km 3 (TJ.) 

Neutrino oscillations can significantly affect the SA^ 
flux, and must be taken into account in the calcula- 
tion of the corresponding lepton event rates in detec- 
tors at the Earth. Earlier studies have stressed, e.g., 
the relevance of — > v T vacuum oscillations @, which 
might produce a few r events per year in a km 3 detec- 
tor (for E > 100 GeV), comparable to the r event rate 
in OPERA, the CERN-to-Gran Sasso long-baseline ex- 
periment [1 61 ] - Although very challenging from an ex- 
perimental viewpoint, the identification of SAv- induced 
r events in future v telescopes might thus supplement 
laboratory tests of — > v T appearance with their astro- 
physical counterpart. Given the potential importance of 
the SAi/ signal detection, we think useful to revisit the 
impact of flavor transitions (in vacuum and in matter) 
on solar atmosphere neutrinos, in the light of the most 
recent advances in v oscillation physics. In particular, 
taking into account updated values for the mass-mixing 
v parameters, we investigate the three-flavor matter ef- 
fects associated to SAz; propagation in the Sun. 

The plan of our work is as follows. In Sec. [Tl] we de- 
scribe the input SAz^ fluxes at the Sun, basically taken 
from Refs. In Sec.|TTT]we study SA^ oscillations in 

vacuum and matter within the standard 3i> mixing frame- 
work, and then calculate their flavor transition probabil- 
ities at Earth by solving numerically the flavor evolution 
equations. In Sec. IIVI we evaluate the impact of flavor 
oscillations on the observable events rates in a km 3 neu- 
trino telescope. We summarize our results in Section fVl 

In a nutshell, we find that interesting, genuine three- 
flavor transitions in matter can take place during the 
SAv propagation, but also that their impact is signifi- 
cantly suppressed by a "conspiracy" of various effects, 
including spatial and energy smearing, and v/V indis- 
tinguishability. In practice, the final oscillation effects 
can be approximated by (phase-averaged) vacuum flavor 
transitions, dominated by the mixing angle #23- 



II. SOLAR ATMOSPHERE NEUTRINO FLUXES 

Cosmic rays hitting the solar atmosphere produce sec- 
ondary particles (mainly it's, plus /x's and ICs), which 
then produce both electron and muon neutrinos and anti- 
neutrinos 1 through their weak decay channels. The most 
recent calculation of the SAv flux has been performed 
in Q by using the Lund Monte Carlo programs PYTHIA 
and JETSET for high-energy particle interactions p7| . 

The fluxes at the source are approximately in the ratio 

cj> e : ^ : <t> T ~ 1 : 2 : , (1) 

as expected by pion production dominance. Moreover, 
an almost equal number of v and V is expected Q , 

v.V~ 1 : 1 . (2) 

After production, SAz/'s are affected by the solar 
medium through both absorption on nucleons and oscil- 
lations in matter. Absorption effects, which play a role 
only at energies E > 100 GeV, depend on the impact 
parameter b of the neutrino trajectory, 

< b < R Q , (3) 

where 6 = corresponds to the Sun center, and b = 
Rq(— 6.96 x 10 5 km) to the solar disk border. In partic- 
ular, for E ~ 10 2 -10 3 GeV, a flux attenuation of 0(10%) 
or less is expected only for cases with b ~ (which have 
little geometrical weight, since > 50% of the flux comes 
from b > 0.7 Rq), while for E > 10 3 GeV absorption 
effects are sizable at any b For E > 100 GeV, the ap- 
proximate flavor-independence of the inelastic neutrino 
cross sections [HI, [l8| allows to factorize neutrino ab- 
sorption from neutrino oscillation effects (which are, in 
principle, entangled) [l9l |. Therefore, in the energy range 
of interest for our work, absorption effects can be taken 
into account through an overall attenuation function mul- 
tiplying the original v fluxes [|[, while oscillation effects 
can be applied afterwards. 

We take the absorption-corrected, unoscillated fluxes 
of v e + V e and 3A. + V ^ from the parametrization in 
Eq. (15) of Ref. [5|, which is integrated over the solar 
disk and is valid for neutrino energies E > 100 GeV. 
For E < 100 GeV, the only calculation of SAis fluxes we 
are aware of is reported in Q, and depends appreciably 
on cosmic ray transport properties in the hclio-magnctic 
fields. Following [||, we assume 4> ^ E~~* , where 7 is al- 
lowed to vary in the range 1.75 < 7 < 2.45 to parametrize 
helio- magnetic uncertainties. 

In Figure Q] we show the unoscillated solar atmosphere 
v fluxes at the Earth, based on the previous input choices. 



1 In the following, we shall sometimes loosely use the symbol 'V 
to indicate both neutrinos and antincutrinos. When needed, a 
distinction between v and v is made to avoid ambiguities. 
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FIG. 1: Unoscillated fluxes of solar atmosphere neutrinos at 
the Earth in terms of neutrino energy. The fluxes are inte- 
grated over the Sun solid angle and summed over the v and 
V components (basically equal). The dominant E~ A depen- 
dence is factorized out. For E > 100 GeV, the fluxes are taken 
from For E < 100 GeV, we take <jj oc E' 1 as in [f|, with 
1.75 < 7 < 2.45 (hatched areas). The overall normalization 
uncertainty (about 20%) is not shown. 



A typical ~ 20% overall normalization error (not shown) , 
due to primary cosmic ray flux uncertainties, should be 
added to the SAi> fluxes. The rapid decrease of the SA^ 
fluxes with energy is only partially compensated by the 
increase of the neutrino cross sections. Figure[5]shows the 
Vft flux (with an intermediate spectral index 7 = 2.1 
for E < 100 GeV) as well as the relevant charged-current 
interaction cross sections with the nucleons [la |. In the 
range 10 < E < 10 3 GeV, the product (p x a decreases 
by more than three orders of magnitude, and thus the 
neutrino event rates above 1 TeV can be ignored in prac- 
tice. For the same reason, relatively low experimental 
thresholds will be crucial for future SAi> detection. 

A final technical comment is in order. We assume a 
^-independent flavor-ratio at the production, which is 
strictly true in the limit in which both pions and muons 
do not lose energy before decaying. We have verified the 
validity of this approximation for muons (for pions, given 
the similar mass but shorter lifetime, the approximation 
would be valid a fortiori). In particular, by using the 
solar atmosphere model reported in we find that, in 
the "worst" case (b ~ and E ~ 1 TeV), the muon range 
before decay is at least 20 times smaller than the typi- 
cal stopping range for a TeV muon in hydrogen (or he- 
lium) [2(|. We estimate that, for E < 1 TeV, the isotropic 
flavor-ratio limit is accurate to better than a few percent, 
in agreement with the results of Ref. [f| . 
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FIG. 2: Energy dependence of unoscillated SA^ fluxes (dot- 
ted line) , of charged current interaction cross sections (dashed 
line), and of their product <f>a (thick solid line) for solar at- 
mosphere + T>i_ l , in arbitrary units. Curves for v c + T> e (not 
shown) are very similar. 



III. SOLAR ATMOSPHERE NEUTRINO 
OSCILLATIONS 

Solar atmosphere neutrinos are affected by flavor os- 
cillations (both in matter and in vacuum) during their 
propagation. Within the standard 2>v framework, the 
oscillation parameters include two squared mass differ- 
ences (5m 2 , Am 2 ) and three mixing angles (612, $23, O13), 
whose ordered rotations form a unitary mixing matrix U. 
We take their best-fit values and 2a ranges from [2l[ : 



Sm 2 


= 7.92(1 ±0.09) x 10~ 5 eV 2 , 


(4) 


Am 2 


= 2.6(l^)xHT 3 eV 2 , 


(5) 


sin 2 6»i2 


= 0.314 (i±g : li) , 


(6) 


sin 2 6*23 


= 0.45 (l±° il), 


(7) 


sin 2 9 13 


= (0.8± 2 ; 3 8 ) x 10- 2 . 


(8) 



The parameter Am 2 has a physical sign, corresponding to 
the cases of normal hierarchy (+Am 2 ) or inverted hierar- 
chy (—Am 2 ). In the following, unless otherwise specified, 
we assume normal mass hierarchy, and also set a possible 
CP-violating phase <5cp to zero. Concerning the mixing 
parameter sin 2 #13, we shall use representative values be- 
low the current upper limits. 

In matter, the v oscillation dynamics also depends on 

[ME!, 



the v e — interaction energy difference V 

±V(x) =±V2G F N e {x) (+/ - for v/v) , (9) 
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FIG. 3: Neutrino potential V(x) in the Sun for three values 
of the impact parameter b. The v trajectory coordinate x is 
normalized to the solar radius. The potential V is shown both 
in units of eV 2 /GeV (left axis) and in units of inverse solar 
radius (right axis). 



where N e is the electron density, which we take from 
the standard solar model in pij , and x is the neutrino 
trajectory coordinate in the range x <E [— a; max , +x max ], 
with i m „ = (Rq — 6 2 ) 1 / 2 . The effective potential V 
encountered by solar atmosphere neutrinos depends on 
the impact parameter < b < Rq . Figure [3] shows rep- 
resentative profiles of V(x). Notice that the maximum 
of V decreases by almost three orders of magnitude as b 
increases from to 2/3i? Q . 



A. The "vacuum" approximation 

Let us provisionally ignore matter effects by setting 
V(x) = in the Sun ("vacuum" approximation, as stud- 
ied in [H). The oscillated SAv fluxes at the detector D 
are given by 



= (pePe, 



^ePe^i 



Vfj,Pfj,fj, 1 



(10) 

(11) 

(12) 



where P a p = P(y a — > z/g). The two independent vac- 
uum neutrino oscillations wavelengths, given by Al = 
AnE/Sm 2 and X H = 4nE/Am 2 , read 



Re 
Ah 

Rq 



4.51 



0.14 



E 



100 GeV 

E 

100 GeV 



7.9 x 10~ 5 eV 2 

5m 2 
2.6 x 10~ 3 eV 2 



(13) 
(14) 



These wavelenghts arc typically much smaller than the 
Sun-Earth distance (L ~ 215 Rq); indeed, in the worst 
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case (E ~ 1 TeV) one has 



B. Matter effects 



L/X L 
L/X h 



5 , 
150 . 



(15) 
(16) 



The occurrence of many oscillation cycles in vacuum 
(L/Xl : h ^> 1), smeared out by realistic energy resolution 
effects (AE/E > 10% is expected for km 3 -sizcd detectors 
[25|). guarantees the incoherence of the SA^ beam at the 
Earth (also when matter effects in the Sun are included, 
as done in the next Sections). 

The vacuum oscillation probabilities are then given by 
their phase-averaged expressions, 



nvac 



i=l 



\U ai \ 2 \Uip\ 2 



(17) 



If we further take sin 2 #13 ~ (a reasonable approxima- 
tion in the context of vacuum oscillations) , then 



In the reference energy range E £ [10, 1000] GeV, so- 
lar atmosphere neutrinos can actually experience peculiar 
three-flavor oscillation effects in matter which, in princi- 
ple, may spoil the "vacuum" approximation discussed in 
the previous section. Matter effects introduce a charac- 
teristic refraction length X r = 2tt/V, which can be ex- 
pressed as 



A r 

Rq 



1.78 x 10" 



eV 



! /GcV \ 



(30) 



V J 

In terms of A r , the oscillation wavelengths in matter read 

Xl 



X, 



X H = 



(cos 20i 2 T Xl/Xt) 2 + sin 2 26>i 2 
X H 

(cos 26»i3 T X H /X r ) 2 + sin 2 20 13 



(31) 



(32) 



pvac 

pvac 
MM 

pvac 
M T 

pv&c 



1 — 9 s 2 r 2 

r, 2 2 2 
/c 12 s 12 c 23 ' 

4(1 - 2s 2 2 c 2 2 ) + a 
2c 2 - 2 n - 2 - 2 



J 23°23 



(1 



'12 l -12 I 



2r 2 s 2 s 2 

zc 12 6 12*23 ' 



(18) 
(19) 
(20) 
(21) 
(22) 



where Cy = ci 
ratio 4> e : 4>^ : 
detector read 



9ij and Sij = sin 6^ . For an initial flux 
r ~ 1 : 2 : 0, the oscillated fluxes at the 



3 [l + 2c 2 2S 2 2 (2c 2 3 -l)] , 

*[! — 2s 23 c 23 + s 12 c 12 c 23 (l — 2c 23 )] 



J2c 23 s 23 



2 2 

Ci 9S1 nS 



24)] 



or, numerically, 



4? 



1.043(1°;°™)^ 

0.493(1^?)^ 
0.485 (+°;°»)^ 



(23) 
(24) 

(25) 



(26) 
(27) 
(28) 



where the above central values are calculated for the 
best- fit values of s 23 and s 2 2 , while the (strongly asym- 
metric) errors are induced by the 2a uncertainties of 



'2:! 



[see Eq. (J7])]. This was recently noticed to be 
a common feature for several astrophysical neutrino 
sources [26|, [2?], [28| . The errors induced by the current 
uncertainties on s 2 2 and on s 2 3 are instead much smaller 
and can be neglected (see also Sec. Ill D below). For 
the specific case of 2-3 maximal mixing (s 23 = 1/2), one 
recovers a well-known result for the phase-averaged, vac- 
uum neutrino flavor ratio at detection 12911 . 



1:1:1 (max. mixing) 



(29) 



where the — (+) sign refers to neutrinos (antineutrinos) . 

Potentially large matter effects are expected when one 
of the above denominators is small, ("resonance condi- 
tion," see [HI), i.e., roughly when 



A 



L,H 



X r 



(33) 



In normal (inverted) hierarchy, the resonance conditions 
are realized only in the neutrino (antineutrino) channel. 
The fact that Sm 2 <C Am 2 also implies that the condi- 
tions Xl.h ~ A,- are realized at rather different energies. 
The dynamics of the 3^ system can then be often ap- 
proximated in terms of two 2v sub-systems, governed by 
(X L ,0 12 ) and (Ah, #13) respectively. 

At resonance, the oscillation lengths in matter are nu- 
merically given by 



\ res 
X L 

Rq 
\ res 

Rpi 



0.97 



0.68 




7.9 x 10- 

Sm 2 /eV 5 

2.6 x 10" 3 
Am 2 /eV 2 J \ sin 



0.93 

sin 26*12 

10- 2 



'13 



, (34) 
(35) 



the latter equation holding for small 613. Very different 
effects take place, when a resonance wavelength A ICS is 
much smaller, comparable, or much larger than the solar 
radius Rq. For A rcs <C Rq (low-energy, adiabatic limit) 
one expects to recover the vacuum case with small, adi- 
abatic matter corrections. For A ros <~ Rq, the solar den- 
sity profile structure is fully probed, and no particular 
approximation or limit is expected to work (also because 
the profile is non-monotonic, see Fig. [3]); this situation is 
realized at E ~ 20 GeV (for A 1 ^ 8 ) or at E ~ O(10 2 ) GeV 
(for A£ s ), see Eqs. fl3HE5j|. Finally, for A ros > Rq, one 
expects again to recover the vacuum limit but for a differ- 
ent reason: large wavelenghts do not "resolve" the solar 
profile structure (extreme nonadiabatic limit [23| ) . In our 
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reference range E £ [10, 1000] GeV, interesting matter 
effects are thus expected to emerge at different energies. 

For 0i3 = these effects are simplest to describe, since 
they are controlled only by (Al,#i2), and just one res- 
onance condition (Al ~ A r ) is effective. At low energy 
(E < 10 GeV), it is A' L CS -C R@, and neutrinos oscil- 
late many times in the solar matter, the phase infor- 
mation being lost. The potential changes slowly over 
one wavelength, so that the adiabatic approximation (no 
transition among mass eigenstates in matter) can be ap- 
plied. The flavor transition probabilities arc then given 
by P a /3 = J2i \U a i\ 2 \Ui/3\ 2 , where U is the mixing matrix 
in matter at the production point. However, since solar 
atmosphere neutrinos are produced in very low density 
regions, one can take U ~ U at the origin, 2 thus recover- 
ing the vacuum limit P a p ~ J2i \U a i\ \Ui@\ at low energy 
(E < 10 GeV). At higher energies (say, E > 50 GeV), it 
is instead X r £ s ^> Rq: neutrino wavelenghts can not re- 
solve the details of the solar matter density profile, and 
the vacuum limit applies again [23| . In conclusion, for 
#13 = 0, significant matter effects are thus expected to 
emerge only around E ~ few x 10 GeV. 

For 9 13 > 0, one expects additional matter effects to 
occur besides the previous ones, since a second resonance 
condition (Ajj ~ A r ) becomes effective. These effects, 
driven by (A#,#i3), are expected to build up at ener- 
gies somewhat higher than for #13 = [e.g., the relevant 
energy range is around E ~ 100 GeV for s 2 3 ~ 10~ 2 , 
see Eq. ([35)1]. It should be noted that: (1) the two 
resonances, governed by different oscillation parameters, 
have different widths; (2) each is realized twice along 
the density profile (with partial cancellation effects); and 
(3) they may take place at energies not widely separated 
(e.g., 20 and 100 GeV), and thus may be partly entan- 
gled and non-factorizablc. These complications prevent a 
straightforward generalization or construction of analyt- 
ical recipes for calculations, in the spirit of those derived, 
e.g., in the context of solar neutrinos [30[. For this rea- 
son, in the following section we shall adopt numerical 
techniques in order to obtain quantitative results for the 
oscillation probabilities in the solar interior. 

Finally, we mention that, in our reference energy range 
E £ [10, 1000] GeV, possible matter effects in the Earth 
(before detection) are negligible, except per hap s at the 
lowest energies (E ~ 10-20 GeV, see e.g. |3l[), where 
they would be anyway much smaller than hclio-magnctic 
flux systcmatics (see Fig. [I]). For the specific case of 
IceCube, at the South Pole, neutrino trajectories in the 
Earth would also be relatively short. In conclusion, wc 
can safely ignore Earth matter effects in the following. 



2 At typical SAu production depths [O(10 3 ) km], the density of 
the solar atmosphere is of O(10 -6 ) g cm -3 , and the neutrino po- 
tential is negligible for our purposes [V ~ O(10 -10 ) eV 2 /GeV, 
not even shown in Fig. 3] . Notice the difference with the more fa- 
miliar adiabatic limit for solar neutrinos, where the high density 
at production (in the solar core) makes U ^ U . 



C. Neutrino oscillation probabilities 

To evaluate the SAi/ oscillation probabilities P a p, we 
evolve numerically — through Runge-Kutta routines — the 
v and V flavor evolution equations [23| for many different 
values of the impact parameter b, corresponding to dif- 
ferent solar matter profiles (see Fig. [3]). More precisely, 
we calculate first the amplitudes AQ(y a — > vi) from the 
production point to the Sun exit, where the v^s are the 
mass eigenstates in vacuum. Then we obtain the flavor 
oscillation probabilities P a p at the Earth as 

3 

P a p = Y J \ A ^^u l )\ 2 -\U. lP \ 2 . (36) 

i=l 

The above factorization of probabilities [23| from the pro- 
duction point to the Sun exit ((^©l 2 ) and from the Sun 
exit to the detector at the Earth (|J7| 2 ) is guaranteed by 
the suppression of all interference terms, due to the many 
oscillation cycles in vacuum and to the finite energy res- 
olution. 

From the unitarity property of the probabilities, one 
can express all the P a ps in terms of only four indepen- 
dent quantities [12]. In the following, we choose P ee , P eM , 
Pfifi and P^t as independent set. 

In Figure 2] we plot these probabilities as functions 
of neutrino energy E after gaussian energy resolution 
smearing (with AE/E = 10% for defmiteness), for a rep- 
resentative impact parameter 6=1/3 Rq, and for both 
s i3 = e ft panels) and s 2 3 = 10~ 2 (right panels), as- 
suming maximal 2-3 mixing and normal mass hierarchy 
in both cases. The upper and middle panels refer to neu- 
trinos and antineutrinos, respectively. Since solar atmo- 
sphere u's and Vs have approximately equal fluxes, and 
are basically indistinguishable in neutrino telescopes, we 
find it useful to plot the arithmetic average of the v and 
V oscillation probabilities in the lower panels of Fig. [U 

For sin 2 f?i3 = (left panels of Fig. 0|, the vacuum 
limit is reached for E > 50 GeV, as discussed in the pre- 
vious section. The same limit applies at very low energies 
[E < 10 GeV, not shown). In the vacuum limit, the prob- 
abilities P a /3 are given by Eqs. (fT8rt22)) . which numerically 
yield: P ee ~ 0.57, P eAl ~ 0.22, P m = P^ ~ 0.39. Mat- 
ter effects, as expected, emerge around E ~ 20 GeV, 
and are more pronounced in the neutrino channel (where 
the resonance condition Al ~ A r can be realized). No- 
tice that, at all energies, P w = P MT . This is due to the 
fact that, for sin 2 #13 = 0, v e v^ T oscillations reduce 
to a pure two- family case (driven by Sm 2 and 6*12), and 
the probabilities in matter can be expressed in terms of 
P ee = P e e{Sm 2 ,9i2,V) alone [3l|. In particular, after 
phase averaging, one gets [32j 

Pfifj. = s 23 + PeeC23 > (37) 
Pfir = (1 + Pee)s| 3 cl3 , (38) 

and thus P w = P MT for s| 3 = 1/2. In addition, one gets 
Pen = (1 — Pee)c2 3 [12], and thus "opposite" variations 
of P ef i and P ee , as evident in the left panels of Fig. 2) 
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FIG. 4: Oscillation probabilities for solar atmosphere neu- 
trinos as functions of energy, for sin 2 Ova = (left panels) 
and sin 2 &i3 = 10~ 2 (right panels). We assume maximal 2- 
3 mixing, normal hierarchy and fix the impact parameter at 
b — 1/3-Rq. Probabilities are shown for v (upper panels), v 
(middle panels), and their average (y + V)/2 (lower panels). 



For sin 2 #13 = 10~ 2 (right panels of Fig. [4j), additional 
three-flavor mixing effects, due to a second resonance, are 
seen to emerge in the neutrino channel around E ~ 100 
GeV, with broad features entangled with the previous, 
lower resonance range. As compared with the left panels, 
the vacuum limit is reached at higher energies, but with 
almost the same limiting values for P a p. Most of the 
additional matter effects appear in i/ e -related oscillations, 
while Pfifj, and P^ T change little, the previous equality 
= P^t being only slightly perturbed. 

The two lowest panels in Fig. [4] show that matter ef- 
fects, averaged over v and v components, produce only 
mild variations with respect to the vacuum limit, inde- 
pendently of the value of sin 2 #13 . For inverted mass hi- 
erarchy, matter effects would be more prominent for an- 
tineutrino oscillation probabilities, but their average over 
v and v would not be much different from Fig. |4] There- 
fore, even if we could select a specific solar atmosphere 
neutrino impact parameter 6, the experimental vjv in- 
distinguishability would make it hard to find signatures 
of matter effects, despite the fact that they may be po- 
tentially large on some Pap's. 



FIG. 5: Oscillation probabilities for solar atmosphere neutri- 
nos, averaged over the impact parameter b, in the same format 
of Fig. m 



Integration over the impact parameter b is expected to 
strengthen these conclusions, since trajectories with rela- 
tively high impact parameter b (and thus with short path 
in the Sun and small matter effects) have large geometri- 
cal weight. Indeed, given the limited angular resolution 
of neutrino telescopes (AO > 1° 25]), the solar disk an- 
gle (A9 sun ~ 0.26°) cannot be resolved, and the neutrino 
oscillation probabilities must be averaged over b: 



{P a0 (E)) 



1 



2irb db P afi {E,b) 



(39) 



Figure 5 shows the 6-averaged oscillation probabilities, 
in the same format of Fig. 4. It appears that possible 
features associated to matter effects are strongly smeared 
out in the observable v + v channel. The energy depen- 
dence of all probabilities is rather flat and roughly in- 
sensitive to sin 2 0i3 , with only minor deviations from the 
vacuum limit governed by the parameter s 2 3 [Eqs. (|18l 
[22])]. Further considerations about realistic experimental 
conditions (integration above a given energy threshold, 
small expected statistics) can only strengthen such con- 
clusions. 
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D. Variations of the oscillation parameters 

The previous results show that observable SA^ fla- 
vor transition effects are dominated by the parameter 
sin 2 6*23- Its relatively large experimental errors [Eq. ([7])] 
represent, at present, the main uncertainty associated to 
oscillation effects in calculating both absolute and rela- 
tive SAis fluxes. In order to illustrate this point we choose 
a representative quantity which might be of interest in 
the future (should v telescopes be able to discriminate 
electron events), namely, the electron-to-muon neutrino 
flux ratio at the detector. 

Figure 6 shows the <fi® /(j>® ratio (with fluxes summed 
over their v and v components) as a function of 
the neutrino energy. The three oscillation parameters 
(Am 2 , 8m 2 , sin 2 ^) are fixed at their best-fit values, 
and normal hierarchy is assumed. The upper and lower 
sets of curves are obtained for the ±2<r extremal val- 
ues of sin 2 6*23 from Eq. ([7]). Within each set of curves, 
sin 2 6*i3 is also taken at its ±2er extrema from Eq. 
while the unknown phase Sep is set at either or i 
(for sin 2 6*13 ^ 0). The spread of the curves due to the 
sin 2 823 uncertainties amounts to about ±15%, while the 
variations due to (sin 2 #23 > 8) are roughly a factor three 
smaller. The ±2cr errors of the other oscillation param- 
eters (Am 2 , 8m 2 , sin 2 6>i 2 ), as well as the change from 
normal to inverted hierarchy, would produce even smaller 
variations in the above curves (not shown). In practice, 
one can currently evaluate the effects of oscillation pa- 
rameter uncertainties in terms of sin 2 $23 only. In the 
future, when sin 2 #23 will be determined much more ac- 
curately, one might need to account also for the (now 
subleading) errors of (sin 2 #13, <5cp)- In the very far fu- 
ture, one might even hope that high-statistics samples of 
SAv events could turn these effects from "uncertainties" 
to "obscrvablcs." 
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FIG. 6: Flux ratio (f>? as a function of energy, for normal 
hierarchy. The spread of the curves reflects the ±2cr spread 
of sin 2 #23 and, to a lesser extent, of sin 2 #13 and Sep- The 
spread of the other oscillation parameters Am 2 , 8m 2 , and 
sin 2 #12 (here fixed to their best-fit values), as well as a change 
to inverted hierarchy, would produce only minor variations in 
the above curves (not shown). 



E. Summary of oscillation effects 

We summarize this Section as follows. We have shown 
that solar atmosphere neutrinos are, in principle, subject 
to significant matter effects in the Sun, at least in some 
oscillation channels and for relatively central trajectories. 
These effects can be sensitive to #13 and to both squared 
mass differences Sm 2 and Am 2 , as well as to the neutrino 
mass hierarchy and to Sqp- However, by averaging over 
the neutrino and antineutrino channel, as well as over the 
impact parameter, these effects are significantly reduced. 

For approximate estimates and for practical purposes, 
one can simply adopt the vacuum limit in Eqs. ()18H22)) . 
and take into account only the uncertainties of #23 (the 
spread of the other oscillation parameters being currently 
much less relevant). Numerically, one can take: 

Pee 0.57 , (40) 

P efi ~ 0.43 c 2 3 , (41) 
P m ~ 0.574 + 4 , (42) 



P MT ~ 1.57 44 , (43) 
Per 0.43 4 , (44) 

where, currently, 4 e [°- 36 > °- 61 ] at 2(7 [ Ec l- ©]■ This 
amounts, e.g., to approximate the curves in the lower 
panels of Fig. 5 with flat, horizontal lines. Of course, 
when SAi/ will be detected in future experiments, it will 
be worthwhile to make the best use of the data by fully 
including matter effects in the oscillation analysis, in or- 
der to avoid unnecessary approximations. For the sake of 
accuracy, we do include matter effects in the calculations 
discussed below. 



IV. FLUXES AND EVENT RATES AT EARTH 

Given the unoscillated input fluxes of Sec. II and the 
neutrino oscillation probabilities discussed in Sec. Ill, one 
can finally obtain the SAv fluxes and event rates observ- 
able at the Earth. Figure [7J shows the oscillated + V ^ 
flux at the Earth, including matter effects (the results 
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Oscillated solar atmosphere v flux at Earth 




Solar atmosphere v event rate 
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FIG. 7: Oscillated flux of solar atmosphere neutrinos 

at the Earth, integrated over the solid angle of the Sun. The 
amplitude of the band above 100 GeV corresponds to the 
current ±2cr uncertainties of the mixing parameter #23 ; below 
100 GeV, the dominant uncertainty is instead related to the 
energy spectral index. 

with no matter effects would be, as already mentioned, 
very similar). The amplitude of the band in Fig. 7 in- 
cludes the 2(T uncertainties in the mixing parameter #23 
from Eq. ([7]). Below ~ 100 GeV, it also includes the dom- 
inant uncertainty related to the energy spectral index 7. 
In particular, the upper curve in Fig. 7 is obtained for 
s 23 = 0-61 and 7 = 2.45, while the lower curve is ob- 
tained for S23 = 0.36 and spectral index 7 = 1.75. The 
fluxes of v e + v e and v T + V T (not shown) would be very 
similar to the flux of + V^, since the flavor ratio at the 
detection is roughly 1:1:1 around maximal mixing. To 
all such fluxes, one should add at least a ~ 20% normal- 
ization error (not shown) due to the uncertainties of the 
primary cosmic ray flux, and to the finite statistics of the 
flux simulations in [||. 

Solar atmosphere neutrinos generate events observable 
in large- volume detectors above a certain energy thresh- 
old Eth- The corresponding event rates (for flavor a) are 
given by 

R a = / dE <P°{E) a a {E) -E- L a {E) A , (45) 
JE th rn P 

where p is the density of the medium (p ~ 1 g/cm 3 in 
water or ice), m p is the proton mass (so that p/m p is 



FIG. 8: v^i and v T event rate in a km 3 detector, as a func- 
tion of the detector threshold energy Eth- The v e rate (not 
shown) is very similar to the one of v^. The width of each 
band reflects the uncertainties on the initial fluxes and on the 
mixing parameters. See the text for details. 

the target nucleon number density), a a is the cross sec- 
tion, L a (E) is the range of the produced lepton (or the 
detector thickness h, whichever the larger), and A is the 
effective area of the detector. The charged-current in- 
teraction cross sections a a are taken from [l^]. For a 
IceCube-like detector we take A = 1 km 2 and h = 1 km. 
Concerning the lepton range, we can take L e = h = L T , 
whereas for 

LJE) = max ( j- In E + , h\ , (46) 
M ' \(3p E th + a/j3 J 1 ; 

where the first term in brackets represents the distance 
travelled by a muon before its energy drops below the 
threshold E t h, the parameters a = 2.5 MeV/(g cm -2 ) 
and (3 = 4.0 x 10~ 6 (g cm -2 ) -1 being adopted from 
More accurate calculations of the event rates would re- 
quire MonteCarlo simulations for specific detector set- 
tings, which are beyond the scope of this work. 

Figure [8] shows the results of our calculation for the 
muon and tau lepton event rates (summing over v and 
77), as functions of the threshold energy Eth, including os- 
cillations in matter. (If the vacuum limit were adopted, 
the curves in Fig. 8 would change by < 3%.) The rate 
of electron events (not shown) would be very similar to 
that of muon events. The width of the bands in Fig. 8 
includes, as in Fig. 7, the uncertainties associated to S23 
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and to 7. On top, one should add a systematic normal- 
ization uncertainty of ~ ±20% (not shown). It appears 
that a relatively low energy threshold (say, at or below 
100 GeV) is crucial to get a handful of events per year in 
km-sized detectors. Note that the rate in Fig. 8 includes 
events from all arrival directions, both above the horizon 
( "downgoing" ) and below the horizon ( "upgoing" ) . 

We conclude this Section with some qualitative re- 
marks about future observations of solar atmosphere neu- 
trino events with rates as low as those calculated in Fig. 8, 
which require careful background rejection. 

The atmospheric neutrino background can be partly 
rejected by directional cuts (in the Sun's direction). Un- 
fortunately, the typical lepton scattering angle and de- 
tector angular resolution are larger than the solar disk 
angle. Therefore, despite the fact that the SA^ signal 
dominates over the atmospheric neutrino noise in the so- 
lar disk radius, the signal-to-background ratio drops to 
< 1 in realistic angular bins In addition, rejection 
of the muon background from cosmic rays showers (still 
huge at 10 2 -10 3 GeV) requires the selection of upgoing 
events, with corresponding reduction of the event rates. 
However, even a few SA^-induced muon events per year 
above a threshold of about 100 GeV (versus a prediction 
twice as large in the absence of 023-driven oscillations) 
might make it possible, in a decade or so, the "solar at- 
mospheric neutrino check" of the large 2-3 mixing estab- 
lished by terrestrial atmospheric neutrinos. 

When a SA^ signal will be established, a rough spectral 
analysis at different energy thresholds might also help to 
check the current expectations for energy spectral index. 
In this respect, our work clearly shows that neutrino os- 
cillations cannot be responsible for dramatic variations 
of the SAv rates with energy; if found, such variations 
should then be ascribed to other factors (e.g., unexpect- 
edly large helio-magnetic effects above 100 GeV). 

The uncertainties on the normalization and energy 
spectrum of the SAv flux can also be reduced through so- 
lar gamma-ray data, since they have common production 
mechanisms (e.g., the 7's from ir° decays). The upper 
limit sets by EGRET on the gamma flux above 100 MeV 
(<i> 7 < 2 x 10~ 7 cm~ 2 s _1 ) J33|] ls onr y a factor ~ 3 above 
the nominal predictions of [3| . The galactic diffuse back- 
ground in the same angular bin is about one order of 
magnitude smaller, and with the performances expected 
by the forthcoming GLAST satellite mission [34| , the so- 
lar gamma signal should be easily detectable. This mea- 
surement may be complemented by the detection of the 
inverse Compton scattering photons produced by cosmic 
ray electrons, which probe the more external layers of so- 
lar atmosphere [35| . All in all, from gamma-ray astron- 
omy we expect a relevant boost in the characterization 
of the solar heliosphere, with consequent sharpening of 
the predictions of SAv fluxes. 

Finally, let us briefly mention a technical point. Since 
the Sun works as a target for cosmic rays (thus generating 
SAz/s), it also shields the Earth from cosmic rays com- 
ing from the same solid angle (thus depicting atmospheric 



z/s in that direction). One might then wonder whether, 
within the solar solid angle, the sought SAz/ signal is ac- 
cidentally compensated by the atmospheric v depletion 
due to the Sun's "shadow," thus escaping detection. We 
argue that this is not the case, for at least two reasons: 
(i) the atmospheric v background depletion is almost one 
order of magnitude smaller than the signal one is look- 
ing for; (it) interplanetary magnetic fields displace (and 
smear) the Sun shadow. The exact displacement is dif- 
ficult to predict, but the one measured with muons by 
MACRO [36| amounts to about 0.6° for median energies 
of the proton primaries of ~ 22 TeV. For the lower ener- 
gies we are focusing on, we expect larger displacements, 
which may allow to separate the Sun shadow in the at- 
mospheric v background from the position of the SA^ 
excess. Finally, the position of this dip (and the similar 
one due to the Moon shadow) will be tagged more accu- 
rately in neutrino telescopes through downgoing muons, 
and thus the deficit should be easily corrected for. 



V. CONCLUSIONS AND PROSPECTS 

During the last forty years, the Sun has been the most 
important source of astrophysical neutrinos, which have 
provided us with crucial evidence of neutrino oscillations. 
Neutrino astrophysics is now entering a new era, where 
next generation neutrino telescopes will be capable to 
detect high-energy (extra)galactic astrophysical neutrino 
fluxes from new sources. These z/s would be useful both 



to probe the properties of the sources |37[ and to test for 
new physics beyond the standard three-neutrino mixing 
scenario [38j | . 

The same telescopes might be able to "see" the Sun in 
a new way, namely, as a source of high energy neutrinos, 
produced by cosmic ray interactions in the solar atmo- 
sphere (SAz>). Even if only a handful of events per year 
may be detectable in a km 3 telescope, their theoretical 
fluxes are relatively under control, and thus a measure- 
ment might provide a high-energy, astrophysical probe 
of neutrino flavor mixing. Moreover, the SAz^ flux rep- 
resents an unavoidable background in searching possible 
neutrino fluxes induced by WIMP annhihilation in the 
Sun core. For these reasons, it appears worthwhile to 
characterize the SAz/ flux detectable at the Earth, after 
neutrino oscillations have occurred in the solar matter 
and in vacuum. 

In our work we have studied solar atmosphere neu- 
trino oscillations in detail, in the light of the most recent 
determinations of the neutrino mass mixing parameters, 
and paying particular attention to matter effects (gener- 
ally neglected a priori in this context). We have found 
that, in the energy range E £ [10, 1000] GeV relevant for 
detection, SAv oscillations in matter have in principle a 
rich phenomenology, and can generate effects which de- 
pend on all 3z^ mass-mixing parameters. Unfortunately, 
these effects arc smeared out to a large extent by sum- 
ming over the neutrino and antincutrino channel, and 
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by integrating over the solar disk angle (as well as over 
energy). In practice, oscillations are dominated by the 
"vacuum" effect of the 6*23 mixing angle. Therefore, even 
if matter effects can and do occur, the original flavor ra- 
tio 4> e : (fift : 4> T ~ 1 : 2 : at the source approximately 
yields the well-known ratio (/>,? : <f>® : <f>® ~ 1 : 1 : 1 at 
the detector (for maximal 2-3 mixing), just as in the case 
of phase-averaged vacuum oscillations. 

Further refinements may include improved estimates 
of the input SA^ fluxes and of their uncertainties. In 
recent years, significant progress has been made in de- 
termining the cosmic ray composition and in refining the 
interaction and cascade development models (3^|. More 
detailed solar atmosphere models and future gamma 
ray data might also help to constrain the SAv produc- 
tion mechanism. In particular, observations in gamma 
rays may sharpen our knowledge of helio-magnetic effects 
3, |35|. Earlier estimates of the unoscillated SAv fluxes 
H, UIl, 0] could thus be usefully revisited. This task will 
certainly become mandatory if a SAv signal will be found 



in neutrino telescopes or, even better, if it will become a 
background for signatures of WIMP annihilation in the 
Sun. 
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